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Abstract 

In this paper, we shall be concerned with geometric functionals and 
excursion probabilities for some nonlinear transforms evaluated on Fourier 
components of spherical random fields. In particular, we consider both 
random spherical harmonics and their smoothed averages, which can be 
viewed as random wavelet coefficients in the continuous case. For such 
fields, we consider smoothed polynomial transforms; we focus on the ge- 
ometry of their excursion sets, and we study their asymptotic behaviour, 
in the high-frequency sense. We put particular emphasis on the analysis 
of Euler-Poincare characteristics, which can be exploited to derive ex- 
tremely accurate estimates for excursion probabilities. The present anal- 
ysis is motivated by the investigation of asymmetries and anisotropies in 
Cosmological data. 
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1 Introduction 



1.1 Motivations and General Framework 

In this paper, we shall be concerned with geometric functionals and excursion 
probabilities for some nonlinear transforms evaluated on Fourier components 
of spherical random fields. More precisely, let {T(x), x G S 2 } denote a Gaus- 
sian, zero-mean isotropic spherical random field, i.e. for some probability space 
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(fi,3,P) the application T(x,u) 1 is {3 x B{S 2 )} measurable, B{S 2 ) de- 
noting the Borcl cr-algebra on the sphere. It is well-known that the following 
representation holds, in the mean square sense (see for instance [2B], [35], [2B]) : 

e 

T(x) = ^2a im Y lm (x) = ^Ti(x) , T £ (x) = 2J ae m Yt m (x) . (1) 

ira i m——£ 

where {Yi m (.)} denotes the family of spherical harmonics, and {ai m } the array 
of random spherical harmonic coefficients, which satisfy Ea£ m a> m ' = C(8\ <5™ ; 
here, 6 a is the Kronecker delta function, and the sequence {Ci} represents the 
angular power spectrum of the field. As pointed out in [30 , under isotropy 
the sequence Cg necessarily satisfies ^1,1 ^1 = ET 2 < oo and the random 
field T(x) is mean square continuous. Under the slightly stronger assumption 
5^ £>L (2£ + l)Ci = 0(log~ 2 L), the field can be shown to be a.s. continuous, an 
assumption that we shall exploit heavily below. 

Our attention will be focussed on the Fourier components {Tg(x)}, which 
represent random eigenfunctions of the spherical Laplacian: 

A S 2T e = -£(£ + l)T e ,1=1,2,... 

A lot of recent work has been focussed on the characterization of geometric 
features for {Ti} , under Gaussianity assumptions; for instance [55], |56j studied 
the asymptotic behaviour of the nodal domains, proving an earlier conjecture 
by Berry on the variance of (functionals of) the zero sets of Ti. In an earlier 
contribution, |12j had focussed on the Deject or signed area, i.e. the difference 
between the positive and negative regions; a Central Limit Theorem for these 
statistics and more general nonlinear transforms of Fourier components was 
recently established by [34]. These studies have been motivated, for instance, 
by the analysis of so-called Quantum Chaos (see again [H]), where the behaviour 
of random eigenfunctions is taken as an approximation for the asymptotics in 
deterministic case, under complex boundary conditions. More often, spherical 
eigenfunctions emerge naturally from the analysis of the Fourier components of 
spherical random fields, as in ([1]). In the latter circumstances, several functionals 
of Ti assume a great practical importance: to mention a couple, the squared 
norm of Ti provides an unbiased sample estimate for the angular power spectrum 
Ce, 

E^J Tl(x)dx^ = (2i + l)C e , 

while higher-order power lead to estimates of the so-called polyspectra, which 
have a great importance in the analysis of non-Gaussianity (see e.g. |29j). 

The previous discussion shows that the analysis of nonlinear functionals of 
{T(} may have a great importance for statistical applications, especially in the 
framework of cosmological data analysis. In this number of papers have 

searched for deviations of geometric functionals from the expected behaviour 
under Gaussianity. For instance, the so-called Minkowski functionals have been 
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widely used as tools to probe non-Gaussianity of the field T(x), see [35J and 
the references therein. On the sphere, Minkowski functionals correspond to the 
area, the boundary length and the Euler-Poincare characteristic of excursion 
sets, and up to constants they correspond to the Lipschitz-Killing Curvatures 
we shall consider in this paper, see pQ, p. 144. Many other works have also 
focussed on local deviations from the Gaussianity assumption, mainly exploiting 
the properties of integrated higher order moments (polyspectra) , see [43] , [44] . 

In general, the works aimed at the analysis of local phenomena are often 
based upon wavelets-like constructions, rather than standard Fourier analysis. 
The astrophysical literature on these issues is vast, see for instance [37J , 0SJ and 
the references therein. Indeed, the double localization properties of wavelets 
(in real and harmonic domain) turn out usually to be extremely useful when 
handling real data. 

In this paper, we shall focus on sequence of spherical random fields which 
can be viewed as averaged forms of the spherical eigenfunctions, e.g. 

Pi(x)=Y,K^j)Tt(x) ,i = 1,2,3... 

i 

for b(.) a weight function whose properties we shall discuss immediately. The 
fields {/3j (a;)} can indeed be viewed as a representation of the coefficients from 
a continuous wavelet transform from T(x), at scale j. More precisely, consider 
the kernel 

f Of _|_ i 

^({x,y)) : =^b( w )-^-P i ({x,y)) 

i 

i m=—t 

Assuming that b(.) is smooth (e.g. C°°), compactly supported in [i? _1 ,i3], and 
satisfying the partition of unity property J2j ^ 2 ("Jr) = 1> f° r au ^ > where is 
a fixed "bandwidth" parameter s.t. B > 1. Then ^j((x,y)) can be viewed as a 
continuous version of the needlet transform, which was introduced by Narcowich 
et al. in [38] , and considered from the point of view of statistics and cosmological 
data analysis by many subsequent authors, starting from [8], [31], [42]. In this 
framework, the following localization property is now well-known: for all M G N, 
there exists a constant Cm such that 

{1 + BJd{x,y)\ 

where d(x,y) = arccos((x, y}) is the usual geodesic distance on the sphere. 
Heuristically, the field 
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is then only locally determined, i.e., for £■? large enough its value depends only 
from the behaviour of T(y) in a neighbourhood of x. This is a very important 
property, for instance when dealing with spherical random fields which can only 
be partially observed, the canonical example being provided by the masking 
effect of the Milky Way on Cosmic Microwave Background (CMB) radiation. 

It is hence very natural to produce out of { (3j (x)} nonlinear statistics of great 
practical relevance. To provide a concrete example, a widely disputed theme in 
CMB data analysis concerns the existence of asymmetries in the angular power 
spectrum; it has been indeed often suggested that the angular power {Ci} may 
exhibit different behaviour for different subsets of the sky, at least over some 
multipole range, see for instance [3T], [13]. It is readily seen that 

e{/$(»)} = £*(^)^*. 

i 

which hence suggests a natural "local" estimator for a binned form of the angular 
power spectrum. More precisely, it is natural to consider some form of averaging 
and introduce 

gj- 2 (z) := J ^K((z,x))/3 2 j (x)dx . 

For instance, should we consider the behaviour of the angular power spectrum 
on the northern and southern hemisphere, we might focus on 

9j . 2 (N):= [ K((N,x))$(x)dx , gj - 2 (S) := [ K((S, x))f3 2 {x)dx , 
Js 2 Js 2 

where K((a, .)) :~ I[ .z.]((a, .)) is the indicator function of the hemisphere cen- 
tred on a <G S 2 , and N, S denote respectively the North and South Poles (com- 
pare [21], [43], [10] and the references therein). More generally, we shall be 
concerned with statistics of the form 

g M (z):= [ K{(z,x))H q (P 3 {x))dx , (2) 
Js 2 

where H q {.) is the Hermite polynomial of g-th order; for instance, for q = 3 
these procedures can be exploited to investigate local variation in Gaussian and 
non-Gaussian features (see [44] and below for more discussion and details). 

1.2 Main Result 

The purpose of this paper is to study the asymptotic behaviour for the expected 
value of the Euler characteristic and other geometric functionals for the excur- 
sion regions of sequences of fields such as {gj-, q (-)}, and to exploit these results 
to obtain excursion probabilities in nonGaussian circumstances. Indeed, on one 
hand these geometric functionals are of interest by themselves, as they provide 
the basis for implementing goodncss-of-fit tests (compare [31]); on the other 
hand, they provide the clue for approximations of the excursion probabilities 
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for {gj;q(-)} , by means of some weak convergence results we shall establish, 
in combination with some now classical arguments described in detail in the 
monograph pQ. 

It is important to stress that our results are obtained under a setting which is 
quite different from usual. In particular, the asymptotic theory is investigated in 
the high frequency sense, e.g. assuming that a single realization of a spherical 
random held is observed at higher and higher resolution as more and more 
refined experiments are implemented. This is the setting adopted in [55], see 
also [S], [27] i [S3] i EH] for the related framework of fixed-domain asymptotics. 

Because of the nature of high-frequency asymptotics, we cannot expect the 
finite-dimensional distributions of the processes we focus on to converge. This 
will require a more general notion of weak convergence, as developed for in- 
stance by [15], [17]. By means of this, we shall indeed show how to evaluate 
asymptotically valid excursion probabilities, which provide a natural solution 
for hypothesis testing problems. Indeed, the main result of the paper, Theorem 
1211 provides a very explicit bound for the excursion probabilities of nonGaussian 
fields such as ^j, e.g. 



where gj- q (x) has been normalized to have unit variance, </>(.), "£(.) denote stan- 
dard Gaussian density and distribution function, a > 1 is some constant and the 
parameters A J;9 have analytic expressions in terms of generalized convolutions 
of angular power spectra, see (|2"2"]) . (fT5]) . 

1.3 Plan of the paper 

The plan of the paper is as follows: in Section 2 we review some background 
results on random fields and geometry, mainly referring to the now classical 
monograph pQ. Section 3 specializes these results to spherical random fields, 
for which some background theory is also provided, and provides some simple 
evaluations for Lipschitz-Killing curvatures related to excursion sets for har- 
monic components of such fields. More interesting Gaussian subordinated fields 
are considered in Section 4, where some detailed computations for covariances 
in general Gaussian subordinated circumstances are also provided. Section 5 
provides the main convergence results, i.e. shows how the distribution of these 
random elements are asymptotically proximal (in the sense of [15]) to those of 
a Gaussian sequence with the same covariances. This result is then exploited 
in Section 6, to provide the proof of ([3]). A number of possible applications on 
real cosmological data sets are discussed throughout the paper. 

2 Background: random fields and geometry 

This section is devoted to recall basic integral geometric concepts, to state 
the Gaussian kinematic fundamental formula, and to discuss its application in 



limsup Pr < sup gj :q (x) > u 





(3) 
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evaluating the excursion probabilities. This theory has been developed in a 
series of fundamental papers by R.J. Adler, J.E. Taylor and coauthors (see [48] . 
[50] . [39] , [3]), and it is summarized in the monographs Q], [2] which are our 
main references in this Section (see also [BJ, [7] for a different approach, and 
[5"T] . [T5] . [3] for some very recent developments in this area). 

2.1 Lipschitz Killing curvatures and Gaussian Minkowski 
functionals 

There are a number of ways to define Lipschitz-Killing curvatures, but perhaps 
the easiest is via the so-called tube formulae, which, in its original form is due 
to Hotelling [22] and Weyl [54]. To state the tube formula, let M be an Tri- 
dimensional smooth subset of R™ such that dM is a C 2 manifold endowed with 
the canonical Riemannian structure on R™. The tube of radius p around M is 
defined as 

Tube(M, p) = {x £ R" : d(x, M) < p} , (4) 

where, 

d{x,M) = inf \\x-y\\. (5) 

Then according to Weyl's tube formula (see [Tj), the Lebesgue volume of this 
constructed tube, for small enough p, is given by 

tn 

A„(Tube(M,p)) = 53p n -Wj£j(M) , (6) 

where uij is the volume of the j-dimensional unit ball and Cj(M) is the j th - 
Lipschitz-Killing curvature (LKC) of M. A little more analysis shows that 
£ m (M) = W m (M), the ?n-dimensional Hausdorff measure of M, and that 
Cq(M) is the Euler-Poincare characteristic of M. Although the remaining LKCs 
have less transparent interpretations, it is easy to see that they satisfy simple 
scaling relationships, in that Cj(aM) — a? 'Cj(M) for all 1 < j < m, where 
aM = {x £ R™ : x = ay for some y £ M}. Furthermore, despite the fact 
that defining the Cj via ([6]) involves the embedding of M in R ra , the Cj(M) are 
actually intrinsic, and so are independent of the ambient space. 

Apart from their appearance in the tube formula ([6j, there are a number of 
other ways in which to define the LKCs. One such (non-intrinsic) way which 
signifies the dependence of the LKCs on the Riemannian metric is through the 
shape operator. Let M be an m-dimensional C 2 manifold embedded in R n ; then 

C k (M) = K n , m ,k f f TT(Sl m - k ^l Nx M(^)n n - m -i(d^n m -i(dx), 

(7) 

where, K n , m: k = ^ m Lk)\ (27r)'"- fc )/ 2 ' an< ^ S(N X M) denotes a sphere in the normal 
space N X M of M at the point x £ M. 



6 



Closely related to the LKCs are set functionals called the Gaussian Minkowski 
functional (GMFs), which are defined via a Gaussian tube formula. Consider 
the Gaussian measure, j n (dx) = (2Tr)~ n ^ 2 e~^ ^ 2 dx, instead of the standard 
Lebesgue measure in ©; the Gaussian tube formula is then given by 

ln ((M,p)) = J^^Ml"(M) , (8) 

fe>0 

where the coefficients Ai]." (M)'s are the GMFs (for technical details, we refer 
the reader to Q]). We note that these set functionals, like their counterparts 
in ([6]) can be expressed as integrals over the manifold and its normal space (cf. 
&)■ 



2.2 Excursion probabilities and the Gaussian kinematic 
fundamental formula 

A classical problem in stochastic processes is to compute the excursion proba- 
bility or the suprcma probability 



P I sup f{x) > u 

\x£AI 

where, as before, / is a random field defined on the parameter space M. In the 
case when / happens to be a centered Gaussian field with constant variance a 2 
defined on M, a piecewise smooth manifold, then by the arguments set forth in 
Chapter 14 of [I], we have that 



Pisnpf(x)<u}-E{C Q (A u (f;M))} < O exp -— ^ , (9) 



2<j 2 



where Co(A u (f; M)) is, as defined earlier, the Euler-Poincare characteristic of 
the excursion set A u (f;M) = {x 6 M : f(x) > u}, and a > 1 is a constant, 
which depends on the field / and can be determined (see Theorem 14.3.3 of [T|). 

At first sight, from ([9]) it may appear that we may have to deal with a hard 
task, e.g. that of evaluating E{Co(A u (f;M))}. This task, however, is greatly 
simplified due to the Gaussian kinematic fundamental formula (Gaussian-KFF) 
(see Theorems 15.9.4-15.9.5 in [T]), which states that, for a smooth M C 1* 

E(4(A(/,M))) 

= d lT( i+ / ) F (l r + (| r + (l i) + 1} (^M+AMWliM), 

e.g., in the special case of the Euler characteristic (i = 0) 

dim(A/) 

E{4(Mf;M))}= E (^)- j/2 C f j (M)M]([u,< X )), (10) 
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where Cj(M) is the j-th LKC of M with respect to the induced metric gf given 
by 

gf(Y x ,Z x )=E{Yf(x)-Zf(x)}, 

for X X ,Y X e T X M, the tangent space at x E M. The Gaussian kinematic 
fundamental formula will play a crucial role in all the developments to follow in 
the subsequent sections. 

3 Spherical Gaussian fields 

In this Section we shall start from some simple results on the evaluation of 
the expected values of Lipschitz-Killing curvatures for sequences of spherical 
Gaussian processes. These results will be rather straightforward applications 
of the Gaussian kinematic fundamental formula (|10p , and are collected here for 
completeness and as a bridge towards the more complicated case of nonlocal 
transforms of Gaussian subordinated processes, to be considered later. 

Note first that for a unit variance Gaussian field on the sphere / : S 2 — >• 
R, the expected value of the Eulcr-Poincarc characteristic of the excursion set 

E{£ (A„(/,5 2 ))} 

= C f (S 2 )MZ([u, ^)) + (2^)- 1 / 2 £{(5 2 )^7 ([u, oo))+{2k)- 1 C{(S 2 )M1([ui 00)) , 
for 

/>OC 

M2(lu,oc))= <f>(x)dx, M]([u,oo)) = Hj-^u^iu), 



where <£>(•) denotes the density of a real valued standard normal random variable, 
and Hj(u) denotes the Hermite polynomials, 

H,{u) = (-1)3 (^)r 1 ^(«) and H^(u) = 1 - *(«), 

while Cl(S 2 ) are the usual Lipschitz-Killing curvatures, under the induced 
Gaussian metric, i.e. 

4(S 2 ) := ( ~ 2 "C /2 / Tr(R^/ 2 )Vol gf ; 



here, R is the Riemannian curvature tensor and Vol g i is the volume form, under 
the induced Gaussian metric, given by 

gf(X, Y) := E {X f ■ Yf] = XYE(f 2 ) . 

We recall that Co (M) is a topological invariant and does not depend on the 
metric; in particular, Co(S 2 ) = 2. Moreover, because the sphere is an (even-) 
2-dimensional manifold, £{(S 2 ) is identically zero. 
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As mentioned before, we start from some very simple result on the Fourier 
components and wavelets transforms of Gaussian fields, e.g. the expected value 
of the Euler-Poincare characteristic for two forms of harmonic components, 
namely 

1 t 
T t {x) = ^2 aimXi m {x) and (3j(x) = ^b(^j)T e (x) , 

m=-i e 

the first representing a Fourier component at the multipole t, the second a 
field of continuous needlet /wavelet coefficients at scale j. We normalize these 
processes to unit variance by taking 



Ti(x) 



T t {x) 



47T *- 



, and pAx) 



i % (21+1) 



Expressions for Lipschitz-Killing curvatures of excursion sets generated by Gaus- 
sian fields are basically known in the astrophysical community (see |35j and the 
references therein). Without too much attention on rigour, our results below 
could be derived from these expressions, but nevertheless we decided to report 
them for completeness and to provide complete mathematical proofs, based 
upon equation (|TU)) . 



Lemma 1 We have 



-in 



I> 2 (X)(2i + i)qp;(i) 

Proof. Recall first that, in standard spherical coordinates 

Pi((x,y)) = P e (smft x smft y cos(cf) x - <p y ) + cos ft x cos ft y ) 
Some simple algebra then yields 



if 



if 



sin ft x sin ft. 



-Pi((x,y)) 



^(i) . 



and 



if 



■Pt((x,v)) 







x=y 



sin ft x dft y dcj) x 

The geometric meaning of the latter result is that the process is still isotropic 
under the new transformation, whence the derivatives along the two directions 
are still independent. We thus have that 



dct 



s 2 



a 2 



sin fiydcftydiJ 



-r(a;,y) 

x=y 

-r(x,y) 



a- 



sin I & x d<f> x d'&> i 

a 2 

sin sin ^ y d4> x d4> r 



■r(x,y) 



x=y 

■r(a;,y) 



1/2 



sin ftdftd<f> 
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Now recall that P' t {l) = ; whence the claim is established. ■ 

Remark 2 Note that since the random field fjj is an isotropic Gaussian random 
field, the Lip 'schitz- Killing curvatures of S 2 under the metric induced by the field 
ft j are given by 

C^(S 2 ) = \)C l (S 2 ), 

where Ci(S 2 ) is the i-th LKC under the usual Euclidean metric, and Xj is the 
second spectral moment of j5j (cf. J7|^. This result is true for all isotropic and 
unit variance Gaussian random fields. 

The second auxiliary result that we shall need is Theorem 13.2.1 in pQ, 
specialized to spherical random fields with unit variance. 



Lemma 3 flT$) For the Gaussian isotropic field f3j : S 2 

~2 — 

0, E/3j = 1, /3j € C 2 {S 2 ) almost surely, we have that 

dim(S 2 )-i 



, such that E/3 ■ 



3 



{c i (A u j (x),S 2 ))}= J2 



i + t 
I 



where 



so that 



Here 



' i + £ ' 






£ 







i/2 



ojiUJi ' Wi ~ r(i + 1) ' 



Pi 



(u) = (2n)- l / 2 M]([u,oo)) = (2n)-( £ + 1 ^ 2 H^ 1 (u)e^ 2 / 2 , 



Po {u) = (2 7 r)- 1 / 2 v / 2^(l - $(w))e u l 2 e~ u ~l 2 = (1 - $(u)) , 
Fn ' 2ir F2V ; 



d 



sin ddcf) 
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Proof. We start by recalling that 



E {^} = ^ E {^} 

whence 

e{£o(A u (^-(x),5 2 ))} = {1 - *(«)} C (S 2 ) + X 



E^ 2 (jr)g^^(l) 



2\ i \ Ue /• /o2\ 



vw 



--C 2 (S 2 



2{1-$(m)} + A : _ 4tt 



2{1-$(u)} + 4tt 



vW 3 
E^ 2 (^^'(i)' 



Also 

E{£ 1 (A u (^.(x),5 2 ))}=^ 



1 r 



1=0 



l+£ 
I 



X e / 2 Pl (u)C 1+e (S 2 ) 



X^ 2 Pl (u)C 2 (S 2 ) 



E^ 2 (jr)C^P;(l) ] 1/2 1 e- 2 / 2 
J V2T 



7T 

— x An x 
2 



f E^ 2 (jr)^^(l) ) 



1/2 



72 



Finally 



E 



[c 2 {A u {P 3 (x),S 2 ))} = Po (u)£ 2 (S 2 ) = {!-*(«)} 4tt 



Using the differential geometric definition of the Lipschitz-Killing curvatures, it 
is easy to observe that 



2E 



S 2 )} = E {Zen(9A„(Tf (.), S 2 )} 



where len(dA u (Te (.), 5 2 ) is the usual length of the boundary region of the 
excursion set, in the usual Hausdorff sense, which can also be expressed as 
d{dA u (T e {.), S 2 ). In fact, writing dA u (f, S 2 ) = f- 1 {u} , we have 



E\cf*(dA u (T £ (.),S 2 ))\=Y, 



k=0 



k + 1 

k 



(2n)-^ 2 M k ({u})C^ +1 (S 2 ) 



The Euler-Poincare characteristic is zero here because we are dealing with a 
closed curve; hence we have 



E\c^(dA u (T e (.),S 2 ))} = 



^ 2 V2 / 1(1+1) 
4 2tt 



(2n)- l l 2 l2 e -j=\x 2 ' 2 C 2 (S 2 ) 



4-7T = 27T 



1) 



- U 2 /2 



11 



where we have used Aii({u}) = 2 e . Hence 

v L 11 V Z7T 

E{£,(a^l„(f,(.),S 2 ))} =A-" 2 E{4'(aA„(f,(.),S 2 ))} 

. x-^kr {*±i> } e-"/ 2 . * |ffi±i> }" 2 e-"' 2 = E {len(OMW, S 2 )} . 

which for u = fits with well-known results on the expected value of nodal 
lines for random spherical eigcnfunctions (see |56j and the references therein). 
Likewise 

E{len(dA u (U),S>)} = 2, { ]^ ^' 2 . (11) 

■ 

Remark 4 ./Vote i/iai, m </ws setting 



' 2 " 




' 2 " 


= 1 , 


' 2 ' 







2 
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In the case of spherical eigenfunctions, the previous Lemma takes the fol- 
lowing simpler form; the proof is entirely analogous, and hence omitted. 

Corollary 5 For the field {T^Q} , we have that 

E {£ (A u (Tz(.), S 2 ))} = 2 {1 - $(u)} + 1 ^ ^ (27r )3 47r ' (12) 
E{£ 1 (A ll (f,(.),5 2 ))}=.|^±^} 1/2 e -« 2 / 2j 

and 

E [£ 2 (A u (r £ (.),5 2 ))} = 4^ x {1 - §(«)} . 

Remark 6 T/?,e expression can be viewed as a weighted average o/E {Zen(cM u (T£(.), 5* 2 )}, 
wzift weights provided by we := b 2 (-gj)C , £ 2 ^ 1 . TTie sequence we is summable; in 
a heuristic sense, we can argue that it must hence be eventually decreasing; we 
thus expect that, for I large enough 

Heb 2 {4r)C^P'eW „ £(£+!) 
I> 2 (i,)^2i±i « 2 > 

and hence 

E{len(dA u 3 (.),S 2 )} <^E{len(dA u (fe{.),S 2 )} . 

This fits with the heuristic understanding that the sequence of fields \/3j(-)} , 
representing an average, is smoother than the subordinating sequence {Tg(.)} ; 
the latter is thus expected to have rougher excursion sets, and hence greater 
boundary regions. This heuristic argument can be made rigorous imposing some 
regularity conditions on the behaviour of the angular power spectrum {Ce} ■ 
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Much more explicit results can of course be obtained by setting a more 
specific form for the behaviour of the angular power spectrum {Cg} and the 
weighting kernel 6(.); for instance, in a CMB related environment it is natural 
to consider the Sachs-Wolfe angular power spectrum Ci ~ Gl~ a , some G > 0, 
a > 2, (see [15]). 



4 Gaussian subordinated fields 



4.1 Local Transforms of 

For statistical applications, it is often more interesting to consider nonlinear 
transforms of random fields. For instance, in a CMB related environment a 
lot of efforts have been spent to investigate local fluctuations of angular power 
spectra; to this aim, moving averages of squared wavelet / needlet coefficients are 
usually computed, see for instance |43j and the references therein. Our purpose 
here is to derive some rigorous results on the behaviour of these statistics. 
To this aim, let us consider first the simple squared field 



H 2j (x) := HMx)) 



E/3 2 Jx) 



1 



The expected value of Lipschitz-Killing curvatures for the excursion regions of 
such fields is easily derived, indeed by the general Gaussian kinematic formula 
wc have 



ElQ{A u {H 2 ;S')) 



Y,{^)~ k,2 £k (S 2 )M^ '((-oo, -Vu + T) U (V^TT, oo)) 

fc=0 
2 



k=0 



4(1 - $(Vu + T)) + 



2- E^dfr)^ 



c 2 {s 2 )- 



-(u+l)/2 



/2tt 



-2y/u + 1 



Likewise 



E\C^(A U (H 2 ;S 2 ))\ = ]T(2 



fc=0 



k + 1 



^ii^ 2 )^ ((-oo, - V^TT) u (Vu+T, oo)) 



l{ ] [s 2 )m% ((-oo, -Vu + i) u (V^TT, oo)) 

+ (2^)- 1 / 2 |4 J (S 2 LMf ((-oo, -v^TT) U (v^TTT, oo)) 
2 l E ^2(^)2|±1Q ' V2^ 



2tt( 



-(u+l)/2 
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which implies for the Euclidean LKC 

1/2 



E{£ l( A u (ff 2 ;S 2 ))} = 2n{ ^ V ^//,Vi/ ) e " 



and therefore 



E^ 2 ( ^)^g^;(i) t _ ( „ +1)/2 

e {A(0^(g a ; ^ 2 ))} = ^ | ^tiS^c^ } e ~ iU+1)/2 - 

Finally 

e(4 3 (A u (# 2 ;S 2 )) j> = £^(5 2 )X^((-oo,-V^TT)U(V^TT,oo)) 
entailing a Euclidean LKC 



E{£ 2 (A U (H 2 ; S 2 ))} = 4tt x 2(1 - $(\ATTT)) . 

It should be noted that the tail decay for the Euler characteristic and the bound- 
ary length is much slower than in the Gaussian case. This is consistent with 
the elementary fact that polynomial transforms shift angular power spectra at 
higher frequencies, hence yielding a rougher path behaviour. Likewise, for cubic 
transforms we have 

±1, 



and finally 

)\=*x{ ^Sr^S^ ) 2 (1 " 
entailing a Euclidean LKC 



E 



[c 2 {A u {p]{x); S 2 ))) = 4tt(1 - 



Similar results could be easily derived for higher order polynomial trans- 
forms. However, although such findings may be useful for applications, as mo- 
tivated above we believe it is much more important to focus on transforms that 
entail some form of local averaging, as these are likely to be more relevant for 
practitioners. To this issue we devote the rest of this section and a large part 
of the paper. 
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4.2 Nonlocal Transforms of /3j(.) 

We now consider the case of smoothed nonlinear functionals. We are interested, 
for instance, in studying the LKCs for local estimates of the angular power 
spectrum, which as mentioned before have already found many important ap- 
plications in a CMB related framework. To this aim, we introduce, for every 

xeS 2 

9r, g (x):= [ K((x,y))H q 3 (y))dy ; (13) 

throughout the sequel, we shall assume the kernel K : [0, 1] — > R to be rotational 
invariant and smooth, e.g. we assume that the following finite-order expansion 
holds: 

£UL Of i i 

K({x,y)) =Y,-^KWPt({x,y}) , some L K G N . 

t 

Here, as before we write H q {.) for the Hermite polynomials. For q = 1, we just 
get the smoothed Gaussian process 

gj(x) := 9j-,i(x) = / K((x,y))'f3 j {y)dy . 

The practical importance of the analysis of fields such as gj- q {-) can be moti- 
vated as follows. A crucial topic when dealing with cosmological data is the 
analysis of isotropy properties. For instance, in a CMB related framework a 
large amount of work has focussed on the possible existence of asymmetries in 
the behaviour of angular power spectra or bispectra across different hemispheres 
(see for instance [43], [44]). In these papers, powers of wavelet coefficients at 
some frequencies j are averaged over different hemispheres to investigate the 
existence of asymmetries/ anisotropics in the CMB distribution; some evidence 
has been reported, for instance, for power asymmetries with respect to the Milky 
Way plane for frequencies corresponding to angular scales of a few degrees (such 
effects are related in the Cosmological literature to widely debated anomalies 
known as the Cold Spot and the Axis of Evil, see jlOj and the references therein). 
To investigate these anomalies, statistics which can be viewed as discrctized ver- 
sions of sup x( zg2gj. q (x) have been evaluated; their significance is typically tested 
against Monte Carlo simulations, under the null of isotropy. Our results below 
will provide the first rigorous derivation of asymptotic properties in this settings. 

Our first Lemma is an immediate application of spherical Fourier analysis 
techniques. For notational simplicity and without loss of generality, we take 
until the end of this Section E/3|(x) = 1, so we need no longer distinguish 

between /3j and 

Lemma 7 The field gj(x) is zero-mean, finite variance and isotropic, with co- 
variance function 

, , p Of J- i 

Eig^xMxz)} = Y / b 2 ( — )K*(e)^—C e P e ((x 1 ,x 2 )) . 

i 
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Proof. Note first that 



E{g j {x 1 )g j {x 2 )} = e\ / K ((x u yi ))p j(yi)dyi / K({x 2 ,y 2 ))p j (y 2 )dy : 

Us 2 Js 2 

K(( Xl , 2/0 )A' «x 2 , y 2 ) )E ( yi )/? . (y 2 ) } dyrdy. 



S 2 xS 2 



/• / 9/4-1 

= / x(( a;ij2/1 )W(.T 2 , 2 / 2 »^fe 2 (— )-^QP,((?/i,y2» 

./s 2 xS 2 V B 4?F 

Recall the reproducing kernel formula 

4-7T 

Pd(xi,yi))Pd\yi>y2))dyi = ^—^ Pi{{xi,y 2 )) , 
P tl {{x 1 ,y 1 ))P l ,{{y 1 ,y 2 ))dy 1 = 0,4^4, 

s 2 



whence 



/ i f (( 2;i , yi ))A'((x 2 , y2 ))^& 2 (-|- ) ?l±lQP,((yi,y2» 

/* o/ _i_ i of _i_ i 

/ E -i 1 "^^'^)! ^^^ 2 )p, 2 ((x 2j2/2 )) 

Js 2 xs 2 ti n e 2 n 

f 9/4-1 

= I £ £ E ^(W 2 (<* 2 ,, 2 >) 



Js 2 47r 



= E 62 (^) K W^^E^T^) / ^(^2 ) W»ft«x 1 ,y a ))d» a 
£ *" e 2 77 Js2 

. , / 9/4-1 

= E &2 (Jj) K2 w^ i ^(^i>^)), 

as claimed. ■ 

The derivation of analogous results in the case of q > 2 requires more work 
and extra notation. In particular, we shall need the Wigner's 3j coefficients, 
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which are defined by (for mi + rri2 + m 3 = 0, see [32], expression (8.2.1.5)) 

i + 4 - 4)K4 - 4 + 4)K4 - 4 + 4)H 17 



lx 4 4 

mi m 2 m 3 



J := (-l)^ 1+mi V24Tl 



(4 + m 3 )!(4 - m 3 )! 



(4+4+4 + 1)! 

1/2 



E 



(4 + mi)!(4 - m0!(4 + m 2 )!(4 - m 2 )!_ 

(-1)*(4 +4 + mi - *)!(4 - mi + 2)! 



z!(4 + 4 - 4 - 4K4 + m 3 - *)!(4 - 4 - m 3 + «)! 



where the summation runs over all z's such that the factorials are non-negative. 
This expression becomes much neater for mi = m 2 = m 3 = 0, where we have 



4 4 4 




, for 4 + 4 + 4 odd 



(-1) 



*l+*2-*3 



i+*2-4 J )/2]![(^i-4!+4 J )/2]![(-£i+*2+4s)/2]! 

for £1 + 4 + 4 even 



+ l2-l3)!(ll-l2+l 3 )!(-ll+l2+l3)! 

(e 1 +i 2 +i3+x)< 
(14) 



It is occasionally more convenient to focus on Clebsch-Gordan coefficients, which 
are related to the Wigner's by a simple change of normalization, e.g. 



/~i£ 3 ?™3 



(-l)*»-"»8 ( tx g 2 

mi m 2 



V24 



1 



£3 
-m 3 



(15) 



Wigner's 3j coefficients are elements of unitary matrices which intertwine alter- 
native reducible representations of the group of rotations S*0(3), and because 
of this emerge naturally in the evaluation of multiple integrals of spherical har- 
monics, see for instance [25], Section 3.5.2. As a consequence, they also appear 
in the covariances of nonlinear transforms; for q = 2, we have indeed 

Lemma 8 The field gj- 2 (x) is zero-mean, finite variance and isotropic, with 
covariance function 

^{g r A x i)9rA x 2)} = 



4tt ^ x Bi 
Proof. Note first that 

E{ffj;2(a?i)#,. 2 (x a )} =E {^ 2 K{ 



4tt 



£ 4 4 




Pe(( x x,x 2 )) 



{.{x u yx))H 2 {fiAyi))dyx / K((x 2 ,y 2 ))H 2 ^Jy 2 ))dy, 



'S 2 xS 2 



K({x llVl ))K{{x 2 ,y 2 ))^{H 2 {f3 J {y 1 ))H 2 {r3 J (y 2 ))}dy l dy 2 
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= 2^ ^K(( a;i ,y 1 ))K((x 2 ,y 2 ))|^fo 2 (-|-)^^QA((2/i,y2»| d Vl dy 2 

= 2 A , E^7^^i)^ 1 ((^i,yi))E^7 ± ^ K ^)P f2 ((x 2 , 2/2 )) 
Js 2 xs 2 e± n e 2 77 

^4 

Now recall that 

' p ii(( x i,yi)) p i 3 i(yi,y2))P£A(yi,y2})dyi 

s 2 



(24 + 1) (24 + 1) (24 + 1) 

x / E ^imi(yi)^iTOi(^l)^3m 3 (yi)^3m3(y2)^4m4(yi)^4m 4 (y2)^l 

= ( (4^ a/a 

l (24 + l)(2£ 3 + l)(2^4 + l) 

x E m 3 m 4 )(o ) ^i-i ^1)^3(2/2)^ (y 2 ), 

mim3m4 ^ / ^ ^ 

Likewise 

A, ((^2, 2/2»^ 3 m 3 (y2)^4TO4 (2/2)^2/2 



.5- 



4tt 



2^ 2 i- 1 Js 2 



T / E r ^m 2 (y2)y« 2 m 2 (2:2)^3^3 (y2)^ 4 m4 (?/2)rf2/2 

I / C2 

M)(24 + i)y-/ h h U\fl% k h\ Y , x) 



V 2£ 2 + 1 ^ I m 2 m 3 m 4 / \ 

V m 2 v / v 

Using the orthonormality properties of Wigner's 3j coefficients (see again 
Chapter 3.5), we have 



E 

7713772 

whence we get 



h h u\( h £ 3 u \ C 2 Ai 



mi m 3 TO 4 J \ 771 2 m 3 7774 / (2£j + 1) ' 

^{9j;2{x 1 )g j . 2 (x 2 )} = 



£ <?iC 2 v 7 
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as claimed. As a special case, the variance is provided by 



Remark 9 Because the field {<?j ; 2(-)} ^as finite-variance and it is isotropic, it 
admits itself a spectral representation. Indeed, it is a simple computation to 
show that the corresponding angular power spectrum is provided by 

r ■ o^m ST }?( h \m e * ) (2£l + 1)(2£2 ± :) r r ( l ^ ^ V 
C tij , 2 .= 2« {l)^b (-)b (— ) C ei C i2 ( o o o ) . 

^1^2 

(16) 

/or £ = 1, 2, ... T/iis result will have a great relevance for the practical implemen- 
tation of the findings in the next sections. 

4.2.1 Higher-order transforms 

The general case of non linear transforms with q > 3 can be dealt with anal- 
ogous lines, the main difference being the appearance of multiple integrals of 
spherical harmonics of order greater than 3, and hence so-called higher order 
Gaunt integrals and convolutions of Clebsch-Gordan coefficients. For brevity's 
sake, we provide only the basic details; we refer to [29] for a more detailed 
discussion on nonlinear transforms of Gaussian spherical harmonics. Here, we 
simply recall the definition of the multiple Gaunt integral (see [29], Remark 6.30 
and Theorem 6.31), which is given by 



Q(£i,m 1 ;...£ q ,m q ;£,m):= / Y timi (x)...Y lq7nq (x)Yi m (x)da(x) 



where the coefficients Q{£\, vrii\ ...£ q , m q ; £, m) can be expressed as multiple con- 
volution of Wigner/Clebsch-Gordan terms (see [15]), 



G(£i,mi; ...£ q ,m q ;£,m) = (-1)' 



'(2li + l)...(2l g + l) 
(4tt)9- 1 (2£ + 1) 



X Z^i Wi0f20 o Ai0£ 3 0-" o A 9 _20£,0 Ly £imi^2ni2 Ly A lA t 1 ^ 3 m 3 --- Ly A 9 _2M,-2^9 m <i ' 

Al...A g _2 Ml---Mg-2 

Following also [29], eq. (6.40), let us introduce the shorthand notation 

r ,Ai...A,_ 2 ^0 ^AiO r t\i<A f-ilO pin it p\ / r ,\i.\ q ~itO\ 

°^i0^ 2 o...« 5 o - — £i0^20 Ai0^ 3 o--- a 9 _20£,o i u l c i > •••> H > 1 ) •— \^e 1 oe 2 o...£ q o j 

Al...Ag_2 

(17) 

It should be noted that, from the unitary properties of Clebsch-Gordan coeffi- 
cients 

£C(4,...,^)= £ {<o%} 2 -E{ C a°_ 2 om} 2 = - = 1 - 
e Ai...a,_ 2 l 
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Lemma 10 For general q > 3, the field gj-q{x) is zero-mean, finite variance 
and isotropic, with covariance function 



5> 2 (€) c{h,...,i q ,i) 



k=l 



Proof. We have 



s 2 Js- 



K((x,y 1 ))K((x,y 2 ))H q (f3 J (y 1 ))H q ((3 J (y 2 ))dy 1 dy 2 



= g !|^|^^(^yi>)^((^y2>)|E 62 (Jj)^^^((yi,2/2))| d Vl dy 2 

It is convenient here to view Tg(x), /3j(x) as isonormal processes of the form 



C t Pe((x,y))dW(y) 



21+1 

47T 



QP^((x,y))dW(y) , 



where dW(y) denotes a Gaussian white noise measure on the sphere, whence 



ii..jt q 

X 



P ei ({x,y 1 ))...P iq ({x,y q ))dW(y 1 )...dW(y q ) 



S 2 x...xS 2 



and 



It follows easily that 



1/6 « ^...f, ^ i=l 

JS 2 x...xS 2 



E te;g( 2 i)S , i; 9 ( z 2)} = 



?2 ~ 2 4 ^ 



4tt 



y 6 2 (^)...6 2 (^) 



, ^ — Ct t \Pi 1 ((x 1 ,x 2 ))...P iq ((x 1 ,x 2 ))dx 1 dx 2 

\ i=i ^ 7r 
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Now write 



(2£ 1 + 1)...(2£, + 1) 

(47r)9 



Pi 1 ((xi i x 2 )).-Pe q ((xi,X2)) 



E Y fim 1 ( x l)- Y i qmq (xi)Y l > imi {x2)-Yi qmq {x2) 



so that 

(24 + 1). ..(21, + 1) 
(47r)« 



S 2 xS 2 



P^ 1 ((zi,xi))Pf 2 ((z 2 ,a;2))P£ 1 ((a;i,a;2))...Pf g ((a;i,X2))cixic;x2 



E E G(h,mi : ...£ q ,m q ; 4,^)0(4, m x ; ...£ q ,m q ;£ 2 , fJ, 2 ) j ^ ^ ^ i/h^O^^O^) | 



4tt 
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The general case <? > 3 hence yields (see also [53], Theorem 7.5 for a related 
computation) 



E 5 yx) 



9'E K2 W E C(^,...,^,^ 2 (^)... 6 2 ( ^^^ 



4 ,24 + 1 2L + 1 



4~ 



and 



E te; 9 (^)5j;<j(y)} 



t l x ...l q 

as claimed. ■ 



Remark 11 It is immediately checked that the angular power spectrum of gj- q (y) 
is given by (see jl7\) ) 



4tt 



Ci-,j. q -dwTT* 2 ® E C(i 1 ,...,£ q ,£)H 
tt...i„ fc=i 



4tt 



(18) 



.4s a special case, for q = 2 we recover the previous result |7 

4 SJ o,4 ,(24 + i)(24 + i) 



c w = 2. 2 (4E^(;^ 2 (|j) 



^2 



4tt 



-Ce x Ce 2 







-in 



4 u2 ,4 ,(24 + i) (24 + 1) 



(19) 



tit 2 



BJ ' 4tt 



4tt 



because 



C(4,4,£) = {Q% } 2 = (2^ + 1) 



* 4 4 
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5 Weak Convergence 



In this Section, we provide our main convergence results. It must be stressed 
that the convergence we study here is in some sense different from the standard 
theory as presented, for instance, by [TT], but refers instead to the broader 
notion developed by [TS], [3], see also [T7], chapter 11. 
We start first from the following Conditions: 

Condition 12 The angular power spectrum has the form 

C t = G(£)£- a , £=1,2,..., 

where a > 2 and G(.) is such that 

< c < G(.) < d 0l 

< Cr e- r , r = 1,2,...,M £ N . 



d r 



Condition 13 The Kernel K(-) and the field {/?.,■(■)} are such that, for all 
j = 1,2,3, ... 

VarU^K^x^H^^dy^o-^B-* , for all j = 1,2,... 

and there exist positive constants c\ , ci such that c\ < crj < ci . 

These assumptions arc mild and it is easy to find many examples such that 
they are fulfilled. Hcuristically, let us consider for instance K((x,y)) = 1, so 
that we are really focussing on weighted averages of the field (y))} . It 

is readily checked that {i? q (/3 -(y))} is a polynomial of order ~ 2 q ^> +1 > and we 
can hence consider the following heuristic argument: we have 



K{{x,y))H q {f5 } {y))dy = / H q {P 3 {y))dy 
s 2 Js 2 

where {£jk>^jk} are a se t of cubature points and weights (see [3S], [5]); indeed, 
because the /3j(.) are band-limited (polynomial) functions, this Riemann sum 
approximations can be constructed to be exact, with weights Xjk of order ~ 
B~ 2j . It is now known that under Condition [T2l it is possible to establish a 
fundamental uncorrelation inequality which will play a crucial role in our proof 
below, see also [8] , [24] . [36] ; indeed, we have that for any M £ N, there exist a 
constant Cm such that 



C M q\ 



{l + Bid(t jkl ,£ jk2 )} qM 
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entailing that the terms H q (/3j(£j k )) can be treated as asymptotically uncorre- 
cted, for large j. Hence 



Var I £ H q (^(C jk ))X 3k \ s £ Far {^(/3 J (^ fe ))} 
[/. x J fee.*,- 

because Efce^f \k — 47r. For instance, for q = 2 we obtain 



\2 



Far <{ j {p){y) - l)dyj = Var | j p){y)dy 



{E^ 2 (5l)(2^+1)Q} 2 {E^ 2 (5t)(2^+1)Q} 2 
{E^ 2 (ir)(2^+1)Q} 2 |i?j( 2 - Q )} 2 



as claimed. 



5.1 Finite-dimensional distributions 

The general technique we shall exploit to establish the Central Limit Theorem is 
based upon sharp bounds on normalized fourth-order cumulants. Note that, in 
view of results from [39], this will actually entail a stronger form of convergence, 
more precisely in total variation norm (see [39]). 

We start by recalling that the field can be expressed in terms of the 

isonormal Gaussian process, e.g. as a stochastic integral 



m = £ b (^(y) = £ b ^f^i^ L mv > z))w{dz) ' 

where W(-A) is a white noise Gaussian measure on the sphere, which satisfies 

EW(A)=0, E{W(A)W(B)}= [ dz, for all A,Be B(S 2 ) . 

J AnB 

It thus follows immediately that the transformed process {H q (f3j(.))} belongs to 
the g-th order Wiener chaos, see [39], [40] for more discussion and detailed defi- 
nitions. Let us now recall the definition of the total variation distance between 
the laws of two random variables X and Z, which is given by 

d T v{X, Z) = sup |Pr(W G A) - Pr(X G A)|. 

AeB(R) 



23 



When Z is a standard Gaussian and X is a zero-mean, unit variance random 
variable which belongs to the g-th order Wiener chaos of a Gaussian measure, 
the following remarkable inequality holds for the total variation distance 



d TV (X 7 Z) < J -t—^-ciim 4 (X) , 

see again [33], [10] for more discussion and a full proof. 

Let us now introduce an isotropic zero-mean Gaussian process fj- q , with 
the same covariance function as that of gj [q . Our next result will establish the 
asymptotic convergence of the finite-dimensional distributions for gj. q and 
In particular, we have: 

Lemma 14 For any fixed vector (x-y, ■■■,x p ) in S 2 , we have that 

9j;q(xi) 9j;q( x p) \ ( /j;<?Ol) U,q( x p) 



LTV 



as j 



V Var (9j;g)' ' V Var (9j.g) I ' \ ^Var{gj. q y ' y/Var(gj- q ) 



= o(l), 



Proof. For notational simplicity, we shall focus on the univariate case; also, 
without loss of generality we normalize j3j(.) to have unit variance. In this case, 
the Nourdin-Peccati inequality ([33], 00]) can be restated as: 



1TV 



, N(0, 1) 1 < . /g-i CMm4 (— ^— ) . (20) 



V Var (9r,q) / V 3<? V Var (9j;q)' 



In view of [201 for the Central Limit Theorem to hold we shall only need to study 
the limiting behaviour of the normalized fourth-order cumulant of gj ;q . Let us 
then consider 

cum 4 {g 3 (x)} 

/ K((x,yi))...K{(x,yi))cumi {H q (/3 j (y 1 )),...,H q (^Jy i ))} dy X -dy 4 



< Ci sup / \K((x,y 1 ))...K((x 7 y 4 ))\\p(y 1 , y2 )\' 1 - r 

r=l,..., q -l J{S2}»4 

x \p(y2,y3)\ r \p(y3,y4)\ q ~ r \p{y4,yi)\ r dyi...dy 4 



<C 2 sup / \p( yil y 2 )\ q r \p{y2,yz)\ r \p{.y3,yi)\ q r \p{yi,yi)\ r dy x ...dy A , 

r=l,...,q-lJ{S' 2 } m 

where 

( ) _ Eib 2 (j J )^C i P e (y 1 ,y 2 ) C M 
P(yUV2) XtW&mCt ~ {l + BH{ yi ,y 2 )} M ' 
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in view of (fT2"j) and the uncorrelation inequality provided by [5], see also [2"4"] . 
[36] . Now standard computations yield 

\p(yi,V2)\ r dy 2 < [ \p{yi,y 2 )\dy 2 < [ — — — ^ -717^2 

< / — w dy 2 + f — w dy 2 < CB-V 

J yr-d{yi,yz)<2B-3 {1 + B->d(yx , y 2 )} J y 2 :d{ Vl ,y 2 )>2B-J {l + B^(yi,2/2)} 

Hence, 

C 2 sup / \p{yi,y2)\ q ~ r \p(y 2 ,y3)\ r \p(y3,y4)\ q ~ r \p(y4,yi)\ r dy^.-dy^ < CB~^ 

r=l,...,?-U{S2}® 4 



and 

[ 9j;g( x ) 



curri4 



{ V Var (9rA x )) 
entailing that for every fixed x £ S 2 , 



0(B~ 2j ) 



d "■(;«' w(0 ' 1) )= o(ir^,, ' 

and hence the univariate Central Limit Theorem, as claimed. The proof in the 
multivariate case (see [41]) is entirely analogous and hence omitted for the sake 
of brevity. ■ 

5.2 Tightness 

We now focus on asymptotic tightness for both sequences {gj- q } and {fj- q } ■ We 
shall exploit the following criterion from [23] : 

Proposition 15 ([23]) Let g,j : AI — > D be a sequence of stochastic processes, 
where M is compact and D is complete and separable. Then g 3 g if gj —tf.d.d. 
g and (tightness condition) 

lim lim sup E( sup \9j(x) — gj(y)\ A 1) = . 

h ~^° j^oc d(x,y)<h 

We are hence able to establish the following 
Lemma 16 For every q £ N, the sequence {gj- q } and {fj- q } are tight. 
Proof. For any Xi,x 2 £ S 2 , we have 



< 



\9rM-9f,M\ < I n \K 3 {(x u y))-K 3 ({x 2l y))\ \H q {(3 3 {y))\dy 

L K d{x 1 ,x 2 ) \H q {Pj(y))\ dy , for all j > 1 . 



5 2 
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where we used the Lipschitz bound 

\ K ji(xi,y)) - K j(( x 2,y))\ < supL K d(x 1 ,x 2 ) , 

j 

which holds because Kj(.) is a finite-order polynomial. Then 

E sup \9j;q(xi) - gj;q(x 2 )\) < sup / L K d(x 1 , x 2 )E I Hq^Pj (y)) \ dy 

\d(xi,x 2 )<h I d(x 1 ,x 2 )<hJ S 2 



< L K h J ^{H q {^{y))Ydy 
= 0(h) , 

whence the result for gj :q follows immediately. Likewise, for any xi,x 2 G S 2 , 
we have 



\ Sup \fj;q(xi) - fj ;q (X2] 
I d(xi ,x 2 )<8 



M 


sup 







y^, a tm(fj;q) {Yt m {xi) - Y em (x 2 )} 
tm 

< ^{E\an m (fj. q )\} I sup \{Y em (xi) - Y £m (x 2 )}\ 

£m [d(xi,x 2 )<5 



Now 



and 



sup \{Y £m (xi) - Y tm (x 2 )}\ < cl 2 5 

d{x\,X2)<& 



^{E|a, m (/ J;9 )|} <^ J{E|a, m (/ J;9 )| 2 } =^(2^+1) v /q(/, ; ,) 

tm tm ' £ 

and because K (.) is compactly supported in harmonic space (and hence, again, 
a finite-order polynomial) 



1/2 



1 £ t (2e + i)c t (f J . g )<o(L K ) 



whence 



E I sup \f m (xi) - f r ,q(x2)\ } < CL 3 K S , 

[d(x u x 2 )<5 J 

for some C > 0, uniformly over j, and thus the result follows. ■ 

5.3 Asymptotic Proximity of Distributions 

Our discussion above shows that the finite-dimensional distributions of the 
non-Gaussian sequence of random fields {gj-, q } converge to those of the Gaus- 
sian sequence {fj- q } ; moreover, both sequences are tight. However, the finite- 
dimensional distributions of neither processes converge to a well-defined limit. 
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In view of this situation, we need a broader notion of convergence than the one 
envisaged in standard treatment such as [TT] ; this extended form of convergence 
is provided by the notion of Asymptotic Proximity, or Merging, of distributions, 
as discussed for instance by [TS], [12], [32], an d others. 

Definition 17 ( Asymptotic Proximity of Distribution f!5f , Ji^f , \ 17p Let g n , f n 

be two sequences of random elements in some metric space (X, p) , possibly de- 
fined on two different probability spaces. We say that the laws of g n , f n are 
asymptotically merging, or asymptotically proximal, (denoted as g n =>■ f n ) if 
and only if as n — > oo 

\Eh(g n )-Eh(f n )\^0 , 
for all continuous and bounded functionals h £ Cb{X,M). 

As discussed by |15j , it is possible to provide a characterization of asymptotic 
proximity that extends in a natural way to standard weak convergence results. 

Theorem 18 Assume that g^ q , fj- q £ Cb{K,S), where K is compact and S is 
complete and separable. Then 9j;q,fj;q are asymptotically proximal if and only 
if they are both tight and their finite- dimensional distribution converge, i.e. for 
all n > 1, x\, ...,x n £ K, we have that 

\^{(9r,M,-,9r, q (x n )) G A} - Pr {(f j;q ( Xl ), f jtq (x n )) G A}\ -> , 
for all A e B(R n ). 

In view of the results provided in the previous subsection, we have hence 
established that 

Theorem 19 As j ' — > oo 

9j\q fj;q > 

i.e. for all h~h: C(5' 2 ,IR) — > R, h continuous and bounded, we have 

\Eh(g Jiq )-Eh(f Jiq )\^0 . 
As a simple application of the asymptotic proximity result, we have 

r S u P9m if sup/ f , i 

I 1 + SUp.g j;9 J 1^ 1 + Sup fj-q J 

It should be noted that asymptotically proximal sequences do not enjoy all the 
same properties as in the standard weak convergence case. For instance, it is 
known that the Portmanteau Lemma does not hold in general, i.e. it is not true 
that, for every Borel set such that Pr {gj £ dA} = Pr {fj £ dA} = 0, we have 

\Pr{g 3 £A}-Pr{f J £A}\^0 . 

As a counterexample, it is enough to consider the sequences fj = —j^ 1 
and gj = However, it is indeed possible to obtain more stringent charac- 
terizations when the subsequences are asymptotically Gaussian. We have the 
following 
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Proposition 20 For every A E B(R), we have that 



Pr <^ sup g j;q (x) -Pr sup f j;q e A 

xes 2 ) Ues 2 



. 



Proof. We shall argue again by contradiction. Assume that there exists a 
subsequence j' n such that for some e > 



Pr \ sup g fn , q (x) e A I - Pr i sup f fn . q G A \ 
UeS 2 ) Ues 2 J 



> £ . 



(21) 



By relative compactness, there exists a subsequence j' 7 [ and a limiting process 
g<x>-. q such that 



Pr < sup g 3 ": q (x) £ i > - Pr < sup g x - q G A 
UeS 2 " J UeS 2 







Likewise, consider {j^"} C {j'n} ; again by relative compactness there exist /oo;g 
such that fj"'- q => foo-,q and hence 



Pr 



sup f j: 

xes 2 



(x) € A } - P r s SUp foo:q £ A 



ies 2 







Note that foo;q, goo;q are isotropic and continuous Gaussian random fields, whence 
the supremum is necessarily a continuous random variable, whence no problems 
with non-zero boundary probabilities arise. Now the finite-dimensional distri- 
butions are a determining class, whence the two Gaussian processes foo;q, goo;q 
must necessarily have the same distribution. Hence 



Pr 



sup fj>"; q (x) e A 

x£S 2 



Pr 



SU P 9jH':q 
x£S 2 



(x) G A 



, 



yielding a contradiction with [21] ■ 

This result immediately suggests two alternative ways to achieve the ultimate 
goal of this paper, e.g. the evaluation of excursion probabilities on the non- 
Gaussian sequence of random fields {gj- q } ■ On one hand, it follows immediately 
that these probabilities may be evaluated by simulations, by simply sampling 
realizations of a Gaussian field with known angular power spectrum; for q = 2, 
for example, fj- q is simply a Gaussian process with angular power spectrum 
given by (|19[) . There exist now very efficient techniques, based on packages such 
as HcalPix ([10]), for the numerical simulation of Gaussian fields with a given 
power spectra; here the only burdensome step can be the numerical evaluation 
of expressions like ()19|) , but this is in any case much faster and simpler than the 
Monte Carlo evaluation of smoothed non-Gaussian fields. Therefore our result 
has an immediate applied relevance. 

One can try, however, to be more ambitious than this, and verify whether 
these excursion probabilities can indeed be evaluated analytically, rather than 
by Gaussian simulations. This is in fact the purpose of the next, and final, 
Section. 
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6 Asymptotics for the excursion probabilities 



The purpose of this final Section is to show how the previous weak convergence 
results allow for very neat characterizations of excursion probabilities, even in 
non-Gaussian circumstances. In particular, if we consider the normalized version 
cjj ;q of the field gj- q , such that gj ]q (x) = [Var(gj ;q (x))) gj ]q (x), then the 
following Theorem can be proved. 

Theorem 21 There exists a constant a > 1, such that 



lim sup 



j- 



Pr 



sup gj- q (x) > u \ - {2(1 - $(u)) + 2u(f>(u)\ j:q } 



< 



exp 



where (see (T8j) ) 



21+1 n p'm 
\ _ 2^,1=1 4?r ^e-j-g^eK 1 ) ,22) 

3'1 x-^L 90-1-1 s-i \ I 



.)■<! 



In order to establish the above Theorem, we shall need to fine tune Theorem 
14.3.3 of PQ to our needs. Let us begin with writing fj. q as a mean zero Gaussian 
random field on S 2 whose covariance function matches with that of gj- q \ note 
that in the previous sections, fj :q denoted the Gaussian random field whose 
mean and covariance matched with that of gj ;q , whereas in this section, the new 
fj.q is basically the normalized version of the old one so as to have unit variance. 
Then, for each xo G S 2 , define 

fr!q{x) = 1 _ ^ Tq ) {fm(%) - P(x,x )f j;q {x ) 

-Co«(/ 3;g (a;),^/ J;(3 (a; ))yar(^/ j;9 (x))^/ j;9 (x) 

d d d \ 

-Cov(fj. q {x), . qa J j;q(xo))Var{ , qa J m ( x )) ■ qa J r,<i( x ) f> 



sm tfacp sm tfu(p sin 

where p(x, xq) = E (fj. q (x)fj- q (xo))- Next define, 

fxj = sup E I sup f*° q (x)) 

and, 

a 2 = sup sup Var(f*° q (x)). 

Xo x^Xq 

Then, under suitable regularity conditions on the kernel K, we have the 
following Proposition, whose proof is deferred to the Appendix. 

Proposition 22 Under the assumption that the kernel K appearing in the def- 
inition of gj-q is of the form 

t 

K{x,y) = Y^k~^Pi{{x,y)), (23) 
i=i 
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for some fixed positive integer £, the field f*° satisfies the following: 

JE(/*»(* 2 ) - f-^)f < K W\X2 - xi\, (24) 

where K(l) depends on q and £, but does not depend on j. 

We are now in a position to provide the following 

Proposition 23 For large enough u, there exists a > 1 such that, uniformly 
over j 

Pr ( sup f m (x) > u\ - EC Q (A u (f m ,S 2 ) < exp(-^-), 

where E£ (A u (f j;q , S 2 ) = 2(1 - + 2u<f)(u)\ j;q . 

Proof. From p. 371 of [1], we know that for u > /it 



Pr { sup fj- q {x) > u 

xeS 2 



Kue 



EC {A u (f m ,S 2 ) 



E|det(-V 2 /^-/i ; ^2) 



1/2 



(25) 



where I2 is the 2x2 identity matrix, K is a constant not depending on j, and detj 
of a matrix is the sum over all the z-minors of the matrix under consideration. 
Our goal is to get a uniform bound for the right hand side of (f2~5j) . Clearly, 

Y^i=o ^ I detj (— V 2 /j :q — fj-.qh) \ is bounded above by a universal constant. To 
get a uniform bound for ^ , we shall resort to the standard techniques on 
estimating suprema of a Gaussian random field using the method of metric 
entropy. En passant, a uniform bound on a 2 is also obtained. 

Applying Theorem 1.4.1 of [1] to the random fields fj- q , and assuming that 
the equation (|24p is satisfied, we conclude that the metric entropy of fj- q can 
be uniformly bounded above, which in turn implies that sup^ < 00. 

Similarly, a uniform lower bound on a 2 is also obtained; for the sake of 
brevity, we shall present its proof in the Appendix. Thus for u > sup.- 



Pr <! sup fj- q {x) > u 
xes 2 



iC (A u (f m ,S 2 ) 



< K ue 



where the K appearing here is different from the earlier one. Finally, note that 
the linear term u can also be ignored by choosing a smaller a in the exponent, 
which completes the proof of the theorem. ■ 

Proof of Theorem 1211 From the results of previous section, we know that 
for any fixed u 



lim 

j-Hx> 



Pr < sup gj- q {x) > u > - Pr <^ sup fj ;q (x) > u 
xes 2 J Ugs 2 



0. 



Combining this with Proposition [23l we get the desired result. 
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7 Appendix 

Here we present the proof of Proposition [22] we shall also obtain a uniform 
lower bound on cr|. For notational simplicity and without loss of generality, we 
take the coefficients in (j23|) to be identically equal to one. 

Writing p(x, y) = cav(fj ;q (x), fj-, q {y)), and , dg x as directional derivatives 
at x in the normalized spherical coordinate directions we have 

cov(^(o:x),^(x 2 )) 

= (i - P ( X0 , iofe - fa, ( p ° T1 ' X2) - p(xo ' Xl)p(xo ' X2) 

-cov(f j . q (jc 1 ),de Xo fj- q (jc ))(xsv(fj. q (x2),d 

- COV (/j;?(^l): 9 0, o /j;?(zo))cOV^^^^ 

-p(x , xi)p{x , x 2 ) + p(x Q , xi)p(x , x 2 )p(x , x ) 
+p(xo,X2)<nv(fj. tg (xi),de xo fj-, g (x ))<nv(fj. tg (x^ 
+p(xo,X2)cav(fj. iq (xi),d4 9 Jj. q (x ))cav(f^ 
-cov(f j;q {x2),de x J j;q (x ))cov(de x J j;q (x2)Jj- q (x 1 ))cov(de x J 

+p(xo,x 1 )cov(f J . q (x2),de x Jj- q (xo))cov(de x Jj- q (x2)J r , q (x ))cov(de x J r , q (x 2 ) 1 de X2 f r , q (x2)) 

+ (var(a exi / i;9 (a;i))) 2 cov(f j - q (x 1 ),de xo fj- q (xo))cov(f j - q (x 2 ),de x Jj- q (x )) 
^v(de x J j - q {xx),d ex J j . q {x 2 )) 

+var(9^ i / j; ,(xi))var(5^ 2 / j;g (a;2))cov(/j; g (xi),9^ o / j;g (a:o)) 
coy(f j;q (x2),de x Jj ;q (xo))coy(de x Jj ;q (x 1 ),d^J j . tq {x 2 )) 

+p(x 07 x 1 )cov(f r:q (x 2 ),d^ x J j . q {x ))var(d^ x J j . q {x2))cov(d^ x J J . q (x2)f r , q {xo)) 
+var(a exi / i;9 (xi))var(a 0x2 / i;g (x 2 ))cov(/ i;9 (xi),a exo / i;9 (a; o )) 

COv(/j ;- q (x 2 ) , 5 0xq fj; q (x ))cov(de xi fj;q{xx),d^Jj. q {x2)) 

+ (varid^Jj^ix!))) cov(f j - q (xi),d^J j;q (xo))cov(f j;q (x 2 ),d^J j;q (xo)) 
cov(9 0xi fj^ixt), d^Jj.^x^fj 

Note that p(xi,x 2 ) can be assumed to have Pi((xi, x 2 )) as the leading poly- 
nomial (uniform over all j). Then, putting x\ = X2 in the above computation, 
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and going through some more (but simple) calculations, one can show that there 
exists a constant M > such that Var(fj°(x)) > M uniformly over all j, which 
in turn, together with the assumption of isotropy, proves that c 2 > M' , for some 
M' > 0. 

Next, to prove the Proposition [22] we begin with 

E (/*°(x 2 ) fZ(*i)) 2 = var(^(xi))+var(^(x 1 ))-2cov(g°(x 1 ) ) ^( a?2 )). 

Wc shall analyze each pair of the terms in the above expression separately 
Let us, for instance consider (together) one of the, seemingly, more involved term 
of the expression which is the last term of the covariance and the corresponding 
term in var(/?° (xi)). At the expense of introducing more notation, let us 
write Ct,4>4, = vax(d^, fj- q (x)) (note that due to isotropy, the variance does not 
depend on the spatial point a;), then the difference between the last term of 
Var(fj.°(xi)) and the last term of Cov{fp{x\), fj.° q (x 2 )), can be written as, for 
all xi, x 2 € (B(xq, e)) c i.e., outside a ball of size e around the point xo, we shall 
have 



(1 - p(x ,xi)) 2 (l - p(x ,x 2 )) 

*( C l,<M> ( COV (f3iq( X l)' d <f>xJr,q( X o))) 0--p(XO,X 2 )) ~ C|^COv(/ J - ;g (xi), fj;q(x )) 
COv(/ j; g(X2),9^ o /j;g(xo))cOv(^ a!i /(x 1 ),9^ 2 /j i .g(x2))(l -p(x ,Xi))) 

(1 - p(x ,xi)) 2 (l - p(x ,x 2 )) 

x (c«;00(l - p(x ,X2))d lj , iiio p(xx,xo) - (1 - p(x , xi))^ p{x 2 , x )d^ d^ 2 p(xi, x 2 fj 

c l<p<p d <f>* P(( x ii x o)) 
(1 - p(x ,xi)) 2 (l - p(x ,x 2 )) 

x ( (?<l>x P( x i> x o) - 9^ o p(x 2 ,x )) Ci-^il - p(x ,x 2 )) 
+ d 4>ro p{x 2l x ) (c w (l - p(x ,x 2 )) - (1 - p(x ,x 1 ))d (j>!ci d^ 2 p(x 1 ,x 2 ) 

Recall that the covariance function p does depend on j, but since we are as- 
suming the kernel K(x, y) to have finite expansion, thus the corresponding Leg- 
endre polynomial expansion of p(xi,x 2 ) can be assumed to have a P^((xi,x 2 )) 
(uniform over j) which is the leading polynomial. Then, taking the modulus of 
the above expression, and considering all x\,x 2 £ (B(xq, e)) c i.e., outside a ball 
of size e around the point xq , we shall have 
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< 



c ld>d> d <t>* n p t(( x i> x o)) 



[1 - P e ((x Q ,x 1 ))]^[l - P i ((x ,x 2 ))} 
x ( ld^ o Pe({xi,xo)) - d^ o Pe({x 2 ,xo))^ C tM [\ - P i ({x ,x 2 ))} 
+ d<f> xo p e({ x 2, x o)) {Q^t 1 - Pi({xo,x 2 ))} - [1 - PidxoiX^d^d^Ptdxi, *2»}) 

{1-P i ((x ,x 1 })) 2 (l-P l ((x ,x 2 )) 

x (J (Pi((xi,x o })(-sin0 Xl sin(0 xl - <j) xo )) - Pi((x2,x o })(-sin0 X2 sin(0 X2 - XO ))) | ■ 
+ |^'((x 2 ,xo))(-sin6> X2 sin(0 a , 2 - ^ XQ ))| • |(C £; ^[1 - P?((x , x 2 ))] 



-[1 - P£((x ,xi))]{P"((xi,x 2 )) sin6> Xl sin 6 X2 sin 2 



<?W) +P/((xi,x 2 ))cos( 



< c\ 



^t;^{^^)[{\ p e{{ x U x o})\ ■ \{smO X2 sin(0 X2 - (f> Xo ) - sin0 Xl sin(0 Xl - 

+|(Pi((x 2 ,x )) -Pi((xi,x )))| ■ |sin0 X2 sin(0 X2 - <£ Xo )|} x eC w 
+M 1 (e,f)C w |P^((xo,xi})-P^((xo,x 2 ))|+M 1 (e,f)|l-P^((x ,xi))| 

x|C«,<W> - P"((xi,x 2 )) sin6» Xl sin6> X2 sin 2 (0 2 , 1 - X2 ) - P^'((xi,x 2 )) cos 
Cl^M(eJ)(eC e ^M 2 (£,e)(\ sm6 X2 \ ■ \ sin(K ~ K) ~ sin (<^i ~ K 



< 



^M(e,£)({|P;((.T 1 ,x ))| • |(sin^ 2 sin(^ 2 - & 

(P^((x 2 ,x )) -Pi((xi,x )))| ■ |sin0 X2 sin(0 X2 - (j> Xo )\} x eQ ;W 

+ M 1 (e,f)|l- 



+M 1 (e,f)C, 

x|Ce,4>4> - P"((xi, x 2 )) sin0 Xl sin# X2 

I sin O X2 - 0xo ) - ^(^x! - <Px ) I 
-M 3 (^e)|x 2 -xi|) 



C 2 ^Af(M)(^ M 2 (£, e 

- sin0 Xl 



+ 1 sin(^ xi 
+M{(e,f)|x 2 

+M('(e,i) x IC^-P;^ 



sin X2 1 
6, 



-<t> XB )\ ■ |sin0 X2 

xi | + M('(e,f) • | sin0 Xl sin 9 X2 \ ■ sin 2 (0 xi - X2 ) 

,x 2 ))cos(0 :Ci - X2 )|J 



< C 2 ^Af (e, £) (eC^M a (e, £)M A (e, I) ■ | sin(«^ 2 - 0, 

+M 4 (e,i) • | sin X2 - sin6 Xl \ + M 3 (e,£)\x 2 - xi\ + M[(e,£)\x 2 - xi\ 
+M'{'{e,l) shi 2 (0 X2 - 9 X1 ) + M'{{e,l) ■ \C tM - P' t {{x x ,x 2 ))\ 

+A/ 1 (ra) (6^)|l-cos(^ 2 -0 Xi 



Now note that is precisely equal to P £ '(l), which can be rewritten as 

P^((xi,xi)). Replacing this in the last part of the above expression, we get the 
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following 



[l-P/((x ,a:i))] 2 [l-P/((a:o > x 2 »] 
x ( {d^Pe^xi^o}) - d^ o Pe((x 2 ,xo))^ Ci-^[\ - Pe((xo, x 2 ))] 
+ d 4>mo Pi{(x 2 ,x a )) - Pi({x ,x 2 ))} - [1 - Pf((a;o,a;i))]90 xi 90 x2 P £ ((xi, 

< Q 2 ^Af(e^)( e ^^M 2 (e^)M 4 (e^) • | sin(^ 2 - <f> x J - sin^ - ^)| 

+M 4 (e,^) • | sin 6*^ - sin 0^1 +M 3 (e,^)|x 2 +M{(e,^)|a; 2 
+M'{'{e,l) sm 2 (6 X2 - 6 Xl ) + M{'(e,£) ■ \P' e ((x 1 ,x l ) - P^((x u x 2 ))\ 

+M[ m \e,t)\l~co^ X2 -$ Xl )\) 

< K(e,£)\ Xl -x 2 \ 

By replicating these set of calculations for each pair of terms in E(fj.°(x 2 ) — 
fj° q (xi)) 2 , we conclude that for every x\,x 2 G B(xo,e), 

Wr^-f-^f^KieJ^-x,]. 

Next we wish to extend this to points inside the set B (xo ■, e)\{xo} , but the 
Lipschitz coefficient K needs to be controlled. Observing that K depends on 
e through the distance of points Xi,x 2 to xq, note that cov(f^° q (xi), (x 2 )) 
grows rapidly as either of X\ or x 2 approach xq, whereas when x\ and x 2 simulta- 
neously approach xq , then the expression assumes the form of an indeterminate 
form, for which one can use the standard l'Hopital's rule to get a precise form 
of the expression. Thus, let us first examine the following 



lim vax(ff.°(x)) 



lim 



(l - p 2 (x ,x) + 2p(x o ,x)d 0:i;a p(xo,x)de ;is p(xo,x)dg !c p(x,x) 



x->x„ (1 - p(x ,x)) 2 

+2p(x ,x)d^ o p(x ,x)d (j , a p(x ,x)d^p(x, x) + {dg xo p(x , x^idj^p^x, x)Y 
p(xo,x)} 2 {dlp(x,x)} 3 " 



Let us do the limit computations for just the first term of the variance 
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expression: 



l-p 2 (x ,x) 

lim ri^r 

x-^xo (1 - p{xo,X)Y 

j. -2p(x ,x)d ea; p(xo,x) 
x^x„ (— 2)(1 - p{xo,x))de x p{x ,x) 

p(xo,x)dg p(xq,x) + (de x p(xo,x)Y 



lira 



(1 - p(x ,x))d^p(x ,x) - (d 8x p(x ,x)) 

p(x ,x)d% p(x ,x) + 3de m p(x ,x)d% p(x ,x) 



x^x (1 - p(x ,x)) d^p(x ,x) - 3d 9x p{x ,x)d^p(x ,x) 

p(x Q ,x)d$ p(x ,x) + d 8:i; p(xo,x)d$ p(x ,x) +3d% p(x ,x)d% p(x , x) + 3d g[c p{x , x)<9| p(x ,x) 

lira - - J ' 

x ^ x ° (1 - p(x ,x)) d^p(x ,x) - 4de x p(x ,x)d$ p(x ,x) - 3 (d^p(x ,x)) 



where we have applied l'Hopital's rule at each step (four times), and we note 
that the final expression is indeed a nontrivial, determinate limit. 

We note that we have assumed p(xo, x) = P(,((xq, x}), and hence the deriva- 
tives above have the following form 



d e:D Pi{(xo,x)) = P' t {-) (cos6 x sin6 Xo cos(<j> x - (j> Xo ) 



sm o x cos ( 



d e x p t(( x o,x)) = P"{-) (cos 6^ sm9 Xo cos((/) x - 4> XQ ) - s,m8 x cos6 Xo ) 

+P'(-) (-sin^ sin 9 Xo cos(<j) x - (j> Xo ) - cos 6^ cos0 xo ) 



dlP e ((x ,x)) 

= P'l' ( ■ ) (cos 9 X sin 8 XQ cos (<p x - cj> Xo ) - sin 9 X cos 9 XQ ) 3 

+2P"{-) (cos9 x sin9x cos(<f) x - <f> xo ) - sin 6^ cos^,,) (- sin^ sin0 Xo cos(cj> x - (f> xo ) - cos9 x cos9 Xo ) 
+P'(-) (-cos 6^ sin 6*^ cos(<j) x - cj) Xo ) + sm9 x cos9 X0 ) 



dlP e ((x ,x)) 

= P^ v \-) (cos9 x sin 6^ cos((/) x - <f) Xo ) - sm8 x cos8 Xo ) 4 

+3P"'(-) (cos 9 X sin 9 Xa cos((f) x - <j) xa ) - sin^ cos 6^) (- sin 6^ sin 6^ cos((j> x - <f> xo ) - cos 6^ cos9 Xa ) 
+2P"'(-) (cos6» x sin 6^, cos((f> x - 4> Xf) ) - sm6 x cos9 Xo ) 2 (- sin 6^ sin 6^ cos(<f> x - 4> xo ) - cos 6^ cos6> Xf) ) 
+2P"(-) (- sin 6^ sin 9 Xo cos(4> x - 4> XQ ) - cos9 x cos9 Xa ) 2 
-2P"(-) (cos 6^ sin 6^ cos(0 x - cf> xo ) - sin9 x cos9 X(1 ) 2 

+P"(-) (cos 6^, sin^ cos((/) x — (f) Xg ) — sin^ cos^) (— cos^ sin^ cos^ — <f> XQ ) + sin9 x cos9 X0 ) 
+P'(-) (sin^a; sin 6^ cos(0 :r — <f> ) + cos9 x cos9 XQ ) 
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Thus we conclude that 

Pi((x ,x)) 
5e x Pf((xo,x)) 
dlP e ((x 0l x)) 
d$ m P e ((x ,x)) 
dlM{xo,x)) 



1 



-P'(l) 


2P £ "(i) + p;(i) 



Subsequently, we shall argue that by continuity, and the fact the field fj° q 
appears to be singular at xq, we conclude that for xi, x-i G B(xq, e) and a small 
enough e, 



sup E 

xi,X2&B(x ,e) 



- fe(*o) 2 = (xi j^ xo , x / - g°,(»o) ! 



The limit on the RHS can again be evaluated by applying l'Hopital's rule, 
and thus, the (uniform) Lipschitz behavior is justified. Thereafter, we note that 

by the isotropy of the underlying field fj ;q , the E fsup xg s2\{ Xo } fj° q ( x )j docs 
not depend on xo, and thus we get a uniform (over j and xq) Lipschitz bound, 
as claimed. 
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